STABILITY OF SYMMETRIC AND NONSYMMETRIC FEM-BEM 
COUPLINGS FOR NONLINEAR ELASTICITY PROBLEMS 
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Abstract. We consider symmetric as well as non-symmetric coupling formulations of PEM 
and BEM in the frame of nonlinear elasticity problems. In particular, the Johnson-Nedelec 
coupling is analyzed. We prove that these coupling formulations are well-posed and allow for 
unique Galerkin solutions if standard discretizations by piecewise polynomials are employed. 
Unlike prior works, our analysis does neither rely on an interior Dirichlet boundary to tackle 
the rigid body motions nor on any assumption on the mesh-size of the discretization used. 



1. Introduction & overview 

The coupling of the finite element method (FEM) and the boundary element method (BEM) 
became very popular when it first appeared in the late seventies of the last century. These 
methods combine the advantages of FEM, which allows to resolve nonlinear problems in 
bounded domains, and BEM, which allows to solve problems with elliptic differential oper- 
ators with constant coefficients in unbounded domains. The two methods are coupled via 
transmission conditions on the coupling boundary 

In 1979, Zienkiewicz and co-workers [ZKB79] introduced a non-symmetric one-equation 
coupling which is based on the first equation of the Calderon system and only relies on 
the simple- layer potential 03 as well as the double-layer potential .ft. In 1980, Johnson & 
Nedelec [JN80] gave a first mathematical proof that this coupling procedure is well-posed 
and stable. This coupling is therefore also referred to as Johnson-Nedelec coupling. Their 
analysis relied on Fredholm theory and the compactness of .ft and was thus restricted to 
smooth coupling boundaries. Based on these works, other coupling methods such as the 
one-equation Bielak-MacCamy coupling and the (quasi-symmetric) Bielak-MacCamy cou- 
pling [BM84] have been proposed. The requirement for smooth boundaries is a severe re- 
striction when dealing with standard FEM or BEM discretizations. Moreover, numerical 
experiments in [CES91] gave empirical evidence that this assumption and hence the com- 
pactness of .ft can be avoided. It took until 2009 when Sayas [Say09] gave a first mathematical 
proof for the stability of the Johnson-Nedelec coupling on polygonal boundaries. 
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In the meantime and because of the lack of satisfying theory, the symmetric coupling 
has been proposed independently by Costabel [Cos88] and Han [Han90]. Relying on the 
symmetric formulation of the exterior Steklov-Poincare operator, [Cos88, Han90] proved 
stability of the symmetric coupling. Early works including [Cos88, CS88, Han90, Ste92] used 
interior Dirichlet boundaries to tackle constant functions for Laplace transmission problems 
resp. rigid body motions for elasticity problems. We also refer to the monograph [GH95] for 
further details. 

To the best of the authors' knowledge, the very first work which avoided the use of an 
additional artificial Dirichlet boundary was [CS95], where a nonlinear Laplace transmission 
problem is considered. In the latter work the authors used the exterior Steklov-Poincare 
operator to reduce the coupling equations to an operator equation with a strongly monotone 
operator. Although their analysis avoids an artificial Dirichlet boundary, their proof of 
ellipticity of the discrete exterior Steklov-Poincare operator, and hence of unique solvability 
of the discrete coupling equations, involved sufficiently small mesh-sizes. Bootstrapping the 
original proof of [CS95], this assumption could recently be removed [AFP 12]. The authors 
of [CFS97] then transferred the ideas of [CS95] to nonlinear elasticity problems in 2D. From 
an implementational point of view, however, the symmetric coupling seems not to be as 
attractive as the one-equation coupling methods, since all four integral operators of the 
Calderon system are involved. 

While Sayas' work [Say09] focused on the linear Yukawa transmission problem as well 
as the Laplace transmission problem, Steinbach [Stell] proved stability for a class of linear 
Laplace transmission problems. He introduced an explicit stabilization for the coupling equa- 
tions so that the stabilized equations turn out to be elliptic. However, the computation of 
the stabilization requires the numerical solution of an additional boundary integral equation 
at every discrete level. Of & Steinbach [OS11] improved the results from [Stell], and also 
gave a sharp condition under which the stabilized problem is elliptic. Based on and inspired 
by the analysis of [Say09, Stell], Aurada et al. [AFF + 12] introduced the idea of implicit 
stabilization. They proved that all (continuous and discrete) coupling equations are equiva- 
lent to associated stabilized formulations, even with the same solution. Since the stabilized 
formulations appear to be elliptic, this proves well-posedness and stability of the original 
coupling formulations, i.e. no explicit stabilization is needed or has to be implemented in 
practice. For the Johnson-Nedelec and Bielak-MacCamy coupling, their analysis covers the 
same problem class as [Stell] and moreover extends it to handle certain nonlinearities. For 
the symmetric coupling, the analysis of [AFF + 12] provides an alternate proof for the results 
of [CS95], but avoids any restriction on the mesh-size. 

In the very recent work [Stel2], Steinbach extended the results from [OS11, Stell] to linear 
elasticity problems. We also refer to [GHS12], where stability of the Johnson-Nedelec, the 
one-equation Bielak-MacCamy, and the (quasi-) symmetric Bielak-MacCamy coupling for a 
Yukawa transmission problem is proven. Moreover, they also show that the Johnson-Nedelec 
coupling applied to elasticity problems with interior Dirichlet boundary is stable for certain 
specific material parameters. 

In our work, we consider (possibly) nonlinear transmission problems in elasticity. As a 
novelty, we introduce a general framework to handle both, the symmetric and non-symmetric 
couplings. We transfer and extend the idea of implicit theoretical stabilization from [AFF + 12] 
to the present setting. This allows us to prove well-posedness of the non-stabilized coupling 
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equations, although they seem to lack ellipticity. The basic idea is the following: We add 
appropriate terms to the right-hand side and left-hand side of the equations and prove that 
this modified (continuous or discrete) problem is equivalent to the original problem, even 
with the same solution. This means that a solution of the modified problem also solves the 
original problem and vice versa. Then, we prove existence and uniqueness of the solution 
of the modified problem and, due to equivalence, we infer that the original problem is well- 
posed. As in [CFS97, GHS12, Stel2], our analysis applies to polygonal resp. polyhedral 
coupling boundaries. From our point of view, the advances over the state of art are fourfold: 

• Unlike [CFS97], we do not have to impose any assumption on the mesh-size h in case 
of the symmetric coupling. 

• Unlike [CS88, GH95, GHS12, Han90, Ste92], we avoid the use of an artificial Dirichlet 
boundary to tackle the rigid body motions. 

• Unlike [Stel2], we prove well-posedness and stability of the original coupling equations 
and thus avoid any explicit stabilization which requires the solution of additional 
boundary integral equations. 

• Unlike [GHS12, Stel2], our analysis for the one-equation couplings also covers certain 
nonlinear material laws, e.g. nonlinear elastic Hencky material laws. 

The remainder of this work is organized as follows: In Section 2, we state the nonlinear 
elasticity transmission problem as well as the precise assumptions on the nonlinearity. Fur- 
thermore, we fix some notation and collect some important properties of linear elasticity 
problems and boundary integral operators, which are used throughout the work. 

Section 3 deals with the symmetric coupling. Here, we introduce the concept of implicit 
stabilization, and prove unique solvability of the coupling equations (Theorem 1). We prove 
that the necessary assumption on the BEM discretization is satisfied, if the BEM ansatz 
space contains the piecewise constants (Theorem 2). 

In Section 4, we apply the ideas worked out in Section 3 to the Johnson-Nedelec coupling. 
Moreover, we incorporate analytical techniques from [OS11, Stel2] to our method and prove 
unique solvability under an additional assumption on the material parameters. 

Finally, the short Section 5 analyzes the one-equation Bielak-MacCamy coupling which 
seems not to be as present as the symmetric resp. Johnson-Nedelec coupling in the litera- 
ture. 

2. Model problem 

Throughout this work, Q C M. d (d = 2, 3) denotes a connected Lipschitz domain with poly- 
hedral boundary T = dfl and complement f2 ext = R d \f2. 

2.1. Notation. We use bold symbols for d- dimensional vectors, e.g. x, and vector valued 
functions u : M. d — > M. d . The components of such objects will be indexed, e.g. u = (ui, u 2 ) T . 
For a set of n G N or a sequence of vector-valued objects we use upper indices for each 
element of the set resp. sequence, i.e. {■u- 7 }™ =1 resp. {u- 7 }!*^. 

Let X C R d be a nonempty, measurable set and let L 2 (X) resp. H l {X), H 1 ^ 2 (X) = 
(i/~ 1 / 2 (X))* denote the usual Lebesgue resp. Sobolev spaces. We define (u , v)x '■= J x uv dx 
for u,v £ L 2 (X). For u G H~ l l 2 {Y),v G H 1 ' 2 ^), the brackets (u , v)r denote the continu- 
ously extended L 2 -scalar product. 
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For vector- valued Lebesgue resp. Sobolev spaces we use bold symbols, i.e. L 2 (X) : = 
[L 2 (X)} d resp. H X {X) := [H 1 (X)} d and so on. Then, we define (u , v) x := f x u ■ v dx 

for u,v G L (X). The product space H := if 1 (ft) x H- 1/2 (T), equipped with the norm 

1 /2 

\\(u,<p)\\ n := (IMItfi(n) + II^IIk- 1 / 2 ^)) for G ^> wil1 be used throughout the 

work. Moreover, let e(u) : cr(v) = Yl'j k=i € jk( u ) (T jk( v ) denote the Frobenius inner product 
for arbitrary tensors e,cr, and define (er(u) , e(v))n ■= f n cr(u) : e{y) dx. The divergence 
div(e(-u)) of a tensor is understood component-wise (div(e(n)))j = J2k=i^ e jk( u )/^Xk for 
j = 1, . . . ,d. Finally, we write ||e(u)||| 2 , n) := (e(u) , e(u)) Q . 



2.2. Linear elasticity. As usual, the linear and symmetric strain tensor e is defined 
component-wise by 

fi\ f \ 1 ( du i , ^ Wfc A 

for all u G -tf and j,k = 1, . . . , d. Together with the Young modulus E > and the 
Poisson ratio v G (0, |), the linear stress tensor cr is defined by 

Ev E 
(2) <j~ k \u) = bjh- — -divu H ejk(u) 

for all n, G and j, fc = 1, . . . , d. To simplify notation, one usually introduces the 

so-called Lame constants 

Ev E 
(3a) A := — and \i :- 



;i + i/)(i-2z/) n 2(i + i/) 

With the identity matrix i" G M dxd , the stress tensor cr then satisfies 

cr(u) = Adiv(u) I + 2fie(u) as well as 
div cr(u) = fiAu + (A + ji) Vdiv(u) in 3D, and 

Ev 

div cr(u) = uAu + ( — + u ) Vdiv(w) in 2D. 

v ' ^ V(l + z/)(l — i/) / v 7 

The kernel of the strain tensor e is given by the space of rigid body motions TZd '■— ker(e) 
{v G i3" 1 (fi) : e(v) = 0} which reads 

(4) K2 :=span{ (j) , (?),("^)} for d = 2 
and 

(5) ^ 3 := span < I I , I 1 I , I I , I x x I , I -x 3 I , I } for d = 3. 




Therefore, it holds cr(v) = for all v G 7^ as well. 
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2.3. Nonlinear transmission problem. As model problem, we consider the following 
nonlinear transmission problem in free space 



(6a) — div2le(fi) = / in Q, 

(6b) -div<x cxt (u cxt ) = mQ ext , 

(6c) u - u cxt —u , on T 

(6d) (2le(u) - <r cxt (u cxt ))n = </> , on T, 

(6e) \u cxt (x)\ = 0(l/\x\) for \x\ ->• oo, 



where n denotes the exterior unit normal vector on T pointing from f2 to f2 ext . The nonlinear 
operator 21 : ffif^ —> Kf^m * s used to describe a (possibly) nonlinear material law in Q. Our 
assumptions on the operator 21 and a more detailed description will be given later on in 
Section 2.5. The stress tensor er ext , which corresponds to the linear elasticity problem in the 
exterior domain, is defined as in (2)-(3) with Lame constants A ext ,yU cxt . For given data / £ 
L 2 (Q),u £ H 1/2 (T), and O £ H~ 1/2 (T), problem (6) admits unique solutions u £ i3" 1 (f2) 
and u ext £ i? 1 1 oc (fi ext ) in 3D. This follows from the equivalence to the symmetric coupling 
and its well-posedness, see Section 3. For the two-dimensional case, the two-dimensional 
compatibility condition 

(7) (f, e J ')n + (0 o , e'') r = j = 1,2 

ensures unique solvability. Here, e? are the standard unit normal vectors in M 2 . We refer 
to [HW08] for further details. 

Remark. The radiation condition (6e) can be generalized to 

(8) u cx Xx) = -G(x)a + r + o(\x\ l - A j for \x\ -> oo, 

with r £ 1Z d , a £ M. d , and G(-) being the Kelvin tensor defined in (10) below. Moreover, 
a = J T (T cxt (u ext )ndT. A solution of (6a)-(6d) with (8) is unique. To see this, we stress that 
the pair (u, u cxt ) solves (6a) -(6d) with (8) if and only if the pair (u, u ext ) = (u — r, u cxt — r) 
solves (6a) -(6d) with 

(9) u x \x) = -G(x)a + 0{l/\x\ 1 - d ) for \x\ -> oo 

and vice versa. Our analysis presented in this work still holds true if we replace (u, u ext ) by 
(u,u cx ) in (6a) -(6d) and the radiation condition (6e) by (9). Note that a = implies the 
compatibility condition (7) in 2D. Therefore, the compatibility condition can be dropped in 
2D for a / 0. In general, the constant a is determined by a = j n fdx + j r (p dT, which 
follows from (6a) and (6d). Furthermore, note that \G(x)\ = 0(l/\x\) for \x\ — > oo and 
d = 3. Hence, (9) coincides with (6e) in 3D. □ 

2.4. Boundary integral operators. The fundamental solution for linear elastostatics is 
given by the Kelvin tensor G(z) £ with 
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for all z G !R d \{0} and j, k = 1, . . . ,d, where G denotes the fundamental solution of the 
Laplacian, i.e. 



;ii) g(z) 




for d = 2, 
for d = 3. 



Throughout this work, 03 denotes the simple-layer potential, .ft the double-layer potential 
with adjoint and 20 denotes the hypersingular integral operator. The boundary integral 
operators formally read for x G V as follows: 

(12) QJ0(aj):= J G{x - y)<t>{y)dT y , 

(13) Siv(x) := J ^G(x - 

(14) 2TJt;(aj) := 

where 7™* denotes the conormal derivative with respect to x defined in (21) below. These 
operators can be extended to continuous linear operators 

(15) 2Jg L(/f- 1/2 (r) ;J Ff 1/2 (r)), 

(16) £G L(/f 1/2 (r);if 1/2 (r)), 

(17) # G L(H~ 1/2 (T); H~ 1/2 (T)), 

(18) QUG L(H 1/2 (T);H- 1/2 (T)). 

We summarize some important properties of these operators. In 3D, the simple-layer poten- 
tial is symmetric and elliptic, i.e. there holds 

(19) (0,QJV)r = (V',^0)r and ||0||^ 1/2(r) < (0 , 5J0) r for all 0, t/> G if~ 1/2 (r). 

Thus, ||0 ||ai := (0) 2J0)J/ 2 defines an equivalent Hilbert norm on H- 1/2 (T). In 2D, elliptic- 
ity can be achieved by an appropriate scaling of the domain Q, see e.g. [Ste08, Section 6.7] 
for further details, and we may thus assume that 03 is elliptic. The hypersingular operator 
is symmetric positive semidefmite, i.e. 

(20) {Wv , w) r = (20w , v) r and (Wv ,v) r >0 for all v, w G H 1/2 (T). 

There holds ker(QU) = ker(| + &) = 1l d , see e.g. [Ste08, Section 6.7]. Throughout this work, 
the boundary integral operators 03, .ft, .ft', and 233 are always understood with respect to the 
exterior Lame constants A cxt ,/i ext . We stress that the natural conormal derivative 7™* is 

(21) jfu := a(u)n on T. 
There holds Betti's first formula, cf. e.g. [Ste08, Section 4.2], 

(22) (a(u) , e(«)) n = (Lu , v) a + ^f(u) , v) r , 
with the linear differential operator Lu = — diver (it). 

2.5. Nonlinear material law and strongly monotone operators. We assume 21 to 
be strongly monotone (23) and Lipschitz continuous (24), i.e. there exist constants c mon > 
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and cup > such that 

(23) c mon ||e(w) - e(v)||| 2(n) < (Sle(u) - 2le(u) , e(tt) - e(v)} n , and 

(24) ||2le(w) - St6(«)|U a( n) < c lip ||e(w) - 6(«)|| £9(n) 

for all u, v G H^fi). In the case 2le(-) = cr(-) there holds, cf. [Ste08, Section 4.2], 

(25) |(«t(«) , e(v)) n \ < d \\e(u) \\ L 2 {n) \\e(v) || £ a (n) , 
and 

(26) e(«)>n>C 2 ||e(tt)||i, (n) 

for all u,v E i? 1 (fi), with constants Ci = 6A + 4/x and C*2 = 2/x. 

An example for a nonlinear material law is the nonlinear elastic Hencky material, obeying 
the Hencky- Von Mises stress-strain relation 

(27) 2te(u) := (K - §£( 7 ( 6 (u))))div(u)I + 2/Z( 7 (e(ti)))e(u) 

with K > being the constant bulk modulus and Lame function 7 (e(w)) := (e(u) — 
-jdiv(it)J) : (e(w) — -jdiv(it)J). Here, ju : M>o — > K.+ denotes a function such that the 
operator from (27) satisfies (23)-(24). Further information on the Hencky material law can 
be found in e.g. [CFS97, CS90, Ste92, Zei88] and the references therein. 

2.6. Discretization. Let l~h denote a regular triangulation of Q and let £\ denote a regular 
triangulation of T. Here, regularity is understood in the sense of Ciarlet. We define the local 
mesh-width function h by h\x '■= diam(X) for X <E % resp. X G ££. Moreover, let K% 
denote the set of nodes of %, and let /C£ denote the set of nodes of ££. We stress that the 
triangulation £^ of the boundary T is, in general, independent of the triangulation l~h- 

Usually, one uses the space V P (S^) := {v G L 2 (T) : v\e is a polynomial of degree < 
p for all E G ££} to approximate functions <p £ H~ l l 2 (T) and the space S q (Th) '■= V q {Th) n 
C(f2) to approximate functions u G with g = p + 1. Here, V q (7~h) '■= {v G 

L 2 (f2) : u|t is a polynomial of degree < g}. In Sections 3-5, we may therefore use the 

space %h '■= %h x 34 = {S q {Th)) x (V p (£^)) to approximate functions (w, </>) G "H : = 
JET 1 ^) x if" 1/2 (r). 

3. Symmetric FEM-BEM coupling 

The symmetric coupling of FEM and BEM has independently been introduced by Costabel 
and Han, see [Cos88, Han90] for example. It relies on the use of all boundary integral 
operators from the Calderon projector. For the derivation of the variational formulation of 
the symmetric coupling, cf. (28), we refer to e.g. [CFS97, CS90, GH95] for nonlinear elasticity 
problems and to e.g. [AFF + 12, CS95, GH95] for nonlinear Laplace problems. It is also shown 
in [CS90] resp. in [CFS97] for the two-dimensional case that the symmetric coupling (28) is 
equivalent to the model problem (6). 

3.1. Variational formulation. The symmetric coupling reads as follows: Find (u,<p) G 
n:=H\Q)x H~ 1/2 {T), such that 

(28a) (2le(u) , e(«)) n + (22Ju , v) T + ((£' - \)4> , v) T = (f , v) n + (<f> + 2Hu , v) r , 
(28b) (t/> , (| - $)u + QJ0) r = (V , (| - #)«o>r 
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holds for all (v, i/>) G %. 

To abbreviate notation, we define the mapping b : H x % — > R and the continuous linear 
functional F G "H* by 

6((u,0),(t>,-0)):=<Sl6(u),6(«))n 

+ (QUit , v) T + ((£' - |)0 , v) T + (V , (§ - £)« + 9J0)r 

and 

(30) V) := (/ , v) n + (0 O + Wu , «) r + (i/> , (§ - £)u }r 

for all (it, 0), (i>, i/j) G "H. Then, the symmetric coupling (28) can also be written as follows: 
Find (it, <j>) G H such that 

(31) b((u,(j>),(v,il>)) = F(v, if)) holds for all (v, if)) G U. 

Note that &(•, •) is nonlinear in it only, but linear in v, if), and </>. If we plug in the functions 
(it, 4>) = (v, i/>) = (r, 0) with r G IZd into (29), we observe 

(32) 6((r,0),(r,0)) = 0. 

Therefore, &(•, •) is not elliptic and unique solvability of (31) cannot be shown directly. In 
the following sections, we introduce an equivalent formulation of (31) which even has the 
same solution. Since this equivalent formulation turns out to be uniquely solvable, also (31) 
admits a unique solution. 

The following two theorems are the main results of this section. With an additional 
assumption on the model parameters c mon , A ext , /i cxt these results also hold true for other 
coupling methods, namely the Johnson-Nedelec coupling, cf. Section 4, and the Bielak- 
MacCamy coupling, cf. Section 5. 

Theorem 1. Let T-L h := X h xy h be a closed subspace ofH and assume thaty c y h nL 2 (r) 
satisfies 

(33) wen d \{0}3£ey (£,r) r ^0. 

Then, the symmetric coupling 

(34) b((u, 0), (v, xj))) = F{v, if)) for all (v,ij>) eU 
as well as its Galerkin formulation 

(35) b((u h , d> h ), (v h , ip h )) = F(v h , iff h ) for all (v h , if> h ) G H h 

admit unique solutions (u,(f>) G % resp. (uh, 4> h ) G Hh- Moreover, there holds the Cea-type 
quasi- optimality 

(36) \\(u,(j>) - (u h ,(f> h )\\ H < Ccea min \\(u,(f>) - (v h ,il> h )\\ H . 

The constant Ccea > depends only on Q, 21, 3^o? an d on the Lame constants A ext ,/i cxt . 

Assumption (33) is clearly satisfied if [Vo : = {'P 1 (£h )) d denotes the space of affine functions 

restricted to ££, since 7Zd C (V 1 (S^)) d and one may thus choose £ = r in (33). However, 

we shall also show that the space y : = (V°{SD) is sufficiently rich to ensure (33). This is 
precisely the second theorem we aim to emphasize and prove. Note that the constant Ccea 
does not depend on the mesh-size h if 3^o Q 34 f° r ai l h. 
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Theorem 2. For d = 2,3, the space y^ := (V°(£^)) d satisfies assumption (33). 

The proof of Theorem 1 resp. Theorem 2 is carried out in Section 3.4 resp. Section 3.5. 

3.2. Implicit theoretical stabilization. To prove Theorem 1, we shall add appropriate 
terms to the linear form b(-, ■), which tackle the rigid body motions in the interior domain 
O. These (purely theoretical) stabilization terms are chosen in such a way that they vanish 
when inserting a (continuous resp. discrete) solution of (34). To be more precise, we will 
use (28b) to stabilize the linear form b(-, ■). 

Proposition 3. Let T-L h — X h x y h be a closed subspace ofH. Let {{$ > j )f =1 } C 34, D G N, 
be a set of linearly independent functions. Define 

D 

(37) 6((u,<£),(«,^)):=K(u,^)>^ 
for all (it, <j>), (v, G %h and 

D 

(38) F(v, := F(v, if>) + ^(tf , (± - Si)u ) (? , (± - &)v + QJ^)r 

3=1 

for all (v,i/j) G %h- Then, there holds the following equivalence: A function (it, 0) G %h 
solves 

(39) b({u,<f>),{v,iJ>)) = F(v,i/>) for all(v,^) E H h 
if and only if it also solves 

(40) b((u,<t>),(v,i>))=F(v^) forall(v^)eU h . 

Proof. Step 1. Assume that (it, <f>) G %h solves (39). Firstly, by inserting the test-function 
(0,^) in (39), we get 

{e , (| - Sl)u + 2J0) r = 6((u, 0), (0, e)) = F(0, = (tf , (§ - &)u ) r 

for all j — 1, ... , D. Secondly, we multiply this equation with (£ J , (■= — &)v + QJi/ , )r and 
infer 

(i j , (| - St)u + QJ0)r(^' , (| - £)« + 2J^)r = , (f - a)«o)r(* j , (§ - + ^)r 

Last, we sum up these terms over all j = 1, . . . , D and add the sum to (39) to see that (it, <fi) 
solves (40). 

Step 2. Assume that (it, (pi) G "Hh solves (40). By choosing (v,ip) = (0,^) as a test- 
function in (40), we infer 

D 

, (i - $L)u + 2J0) r + > (I - &) u + ^>r(f , 9J^)r = &((«, 0), (0, 

= F(0, O = , (| - il)tio)r + E^' > (5 - *)«o>rte J ' , QJtf )r, 

j'=i 
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which is equivalent to 

D 

(41) , (i "•&)(« " «o) + 2J0)r(f , 9J^)r = , {\ ~ &){u - u ) + 2J0) r 

3=1 

for all £ = 1, . . . , D. Next, we define a matrix A E IR^ m D witn entries A jfc := (£ fc , 2J^) r 
and a vector a; e R D with entries x k := (£ fc , (| — .ft) (it — it ) + 23</>) r for all j,k — 1, . . . ,d. 
Then, we can rewrite (41) for all ^ = 1, . . . , D simultaneously as 

(42) Ax = -x. 

Since 23 is elliptic, and {£ J )f = i are linearly independent, the matrix A is positive definite 
and thus only has positive eigenvalues. Therefore, (42) is equivalent to x = 0, which means 

(e , (| - &)u + QJ0)r = ($ j , (§ - £)u > r 
for all j = 1, ... , -D. With these equalities and the definitions of •) and &(•, •), we get 

D 

6((u, 0), («, V)) - &((«, 0), (»> </>)) = > (l " ^ + ^Mtf , (I " + ^)r 

i=i 

= > (I - #) w o)r<^ , (| - *)« + 2JV)r 

i=i 

= F(i>,-0)-F(i;,<0). 

In particular, (40) thus implies (39). This concludes the proof. □ 

3.3. Equivalent norm. To show that the equivalent bilinear form of Proposition 3 is, in 
fact, stabilized and yields a strongly elliptic formulation, we will prove that the employed 
stabilization term provides an equivalent norm on the energy space %. 

Lemma 4. Let gj : % — > K with j = 1, . . . , D denote linear and continuous junctionals such 
that 

D 

(43) 1^,0) | 2 := 53&(r,0)V holds for all r E TZ d \{0} . 

3=1 

Then, the definition 

(44) |||K0)||| 2 := ||e(w)||! 2(n) + (</>, 2J0) r + |^K0)| 2 /or all (u,<f>) E H 
yields an equivalent norm on H, and the norm equivalence constant C nor m > in 

(45) C^ m \\(u,ct>)\\u< \\\(u, 0)||| <C aoim \\(u,<f>)\\ H for all (u, cf>) E H 
depends only on Q, A ext ; /i ext ; and g. 

Proof. Firstly, due to boundedness of g and || II x, 2 (n) ~ IMIh 1 ^) there holds \\\(u, tj>)\\\ < 
|| (it, 4>)\\h- Secondly, we argue by contradiction to prove the converse estimate: Assume that 
there exist functions (it n , <p n ) with \\(u n , (f} n )\\n > ^IIK^n, n )lll f° r a ^ n EN. Define 



\\{u n ,4> n )\\ H ' 
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Then, it follows \\\(v n , t/OII < V 71 an< ^ thus e jk(v n ) — > in L 2 (Vt) for j, = as 
well as i\) n — > in i/~ 1 / 2 (r). By definition of (v n ,if> n ), there holds ||(v n , VOIIft = 1 an d 
we may extract a weakly convergent subsequence with (v ne ,tf) ne ) — ^ {v,ij)) in "H. Next, we 
conclude that i/> n — )■ i/> = in i3"~ 1 ^ 2 (r) and v ne — > v in i 2 (f2), where the latter follows 
from weak convergence v ne ^ v in and the Rellich compactness theorem. Weak 

lower semi-continuity of |||(-,-)lll implies ||| (v, if}) ||| = and thus e(v) = and |g(t;,0)| = 0. 
Due to ker(e) = IZd and (43), it follows that v = 0. Moreover, with Korn's second inequality, 
cf. e.g. [Ste08, Theorem 4.17], we infer 

IK, - v\\ 2 H i { n) ~ ll e KJ ~ c(«)||?,3(n) + IK - v \\l\n) 

and therefore (v ne , ij) n ) — >■ (0, 0) in "H, which contradicts IK^, V'njllw — 1- This concludes 
the proof. □ 

The following proposition provides the equivalent norm used to analyze the symmetric 
coupling as well as the Johnson- Nedelec coupling (see Section 4 below). 

Proposition 5. Let y$ C y n L 2 {Y) be a subspace which satisfies assumption (33) of The- 
orem 1. Let r 1 , . . . , r D with D = dim(7?.d) denote a basis of the rigid body motions and let 
Ho : L 2 {Y) — > be the L 2 -orthogonal projection. Then, £ J := Ho^) for j = 1,...,D are 
linearly independent. Moreover, the functionals gj £ Ti* defined by 

(46) gj (u,<f>):=(e,(l-Si)u + W<j>)r for (u,cf>) eH 
fulfill assumption (43) of Lemma 4- In particular, 

D 

(47) |||(u, 0)||| 2 := \\e(u)\\l HQ) + (<f> , 2J0) r + £ |<tf , (§ - *)« + 2J0) r | 2 

an equivalent norm on H, and the norm equivalence constant C norm > in 

(48) C£L ||(«,^)||« < I («, 0)||| < Cnorm ||(u,0)|| w /or a// (u» £ ft 
depends only on Q, 3^o? -^ cxt ; /i ext . 

Proof. We work out an alternative formulation of (33). With , r)p = (n £ , r)r = 
(£ , n r)r, condition (33) becomes 

Weii d \{o}3£ey (£,n r) r ^o. 

Clearly, this is equivalent to IT (r) ^ for all r £ 7£d\{0}, which yields that the functions 
£ J := no(W), for j = are linearly independent. Therefore, we can reformulate 

condition (33) as 

(49) Vr£^\{0}3j£{l,..., J D} (? , r) r ^ 0. 

The functionals are well-defined, linear, and bounded. To see (43), we stress that due to 
ker(| + £) = ft d , 

g 3 (r, 0) = , (f - £)r} r = , r) r for j = 1, . . . , D and r £ 7^. 

n 



From (49) we infer that there exists j G {1, . . . , D} such that gj(r, 0) ^ 0. Therefore, (43) 
holds for 

D D 

\g(u, 0)| 2 = 0) 2 = , (i - *)« + 5J0)r| 2 - 

This concludes the proof. □ 

3.4. Proof of Theorem 1. As far as existence and uniqueness of solutions is concerned, it 
suffices to consider the Galerkin formulation (35), since this covers the case Hh = Ji as well. 
With assumption (33) and y Q 34 H L 2 (T), Proposition 5 allows to apply Proposition 3. 
Hence, we may equivalently ask for the unique solvability of (40) instead of (39) resp. (35). 
To this end, we define the nonlinear operator 03 : %h H* h by 



h)i 



First, we rewrite equation (40) as an equivalent operator equation: Find (uh, </>/J £ Hh such 
that 

(50) %(u h ,(f> h ) = F mU* h . 

Step 1 {Lipschitz continuity o/QS). Due to the Lipschitz continuity (24) of 21 and the 
boundedness of the boundary integral operators, it clearly follows that *B is also Lipschitz 
continuous. The Lipschitz constant C\\ v > in 

(51) < OipIKu,,,^) - K,^)lk 

for all (it/j, (vfe, ■0 h ) G "H, thus depends only on 21, f2, A cxt , and yU. ext . 

Step 2 [Strong monotonicity o/53). We have to prove that, for all (u^, h ), (u^ ip h ) G "H, 

(52) (»(«/», (f> h ) -*8(v h ,i/> h ), (u h -v h ,(p h -^ h )) > C mon \\(u h - v h ,(j> h - Vh)llw- 
To abbreviate notation, let (tw^, x h ) := (it^ — u^, </) h — Then, we get 

= (Sae(tih) - ae(uft) , e(w h ))a + (2Uu> ft , i^) r + ((# - \)Xh > ™fc)r 



5/ 1' 

2' 



+ ( Xfc , (| - + K Xh )r + 52\(Z j , (| - &)w h + OJx^rl 2 =: I 

Next, we use strong monotonicity (23) of 21 and positive semi-definiteness (20) of Q2J to 
estimate 

D 

I > c mon \\e(w h )\\l 2 ^ + ( Xh , K Xh )r + ' (§ " ^ Wh + *Xh)r\ 2 

D 

>min{c mon ,l} (\\e{w h )\\ 2 LHn) + ( Xh ,V3x h )r + J2\(€ J > d~ &) w h + Wx h )r\ 
= min{c mon ,l} HKw^Xfc)! 2 - 
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Finally, the norm equivalence of Proposition 5 yields strong monotonicity, where C mon = 
min{c mon , l}Cnorm > depends only on 21, fl, A cxt , /x ext and 3V 

Step 3 (Unique solvability and Cea lemma). The main theorem on strongly monotone 
operators, see e.g. [Zei90, Section 25], states that the operator formulation (50) and thus 
the Galerkin formulation (35) admits a unique solution (uh, 4>h) e Hh- F° r Hk — H, we 
see that also the symmetric formulation (34), admits a unique solution (it, 4>) G H. Finally, 
standard theory [Zei90, Section 25] also proves the validity of Cea's lemma (36), where 
Ccea = Ciip/Cmon > depends only on Q, 21, A cxt , /i oxt and y . □ 

Remark. Our analysis unveils that (28b) tackles the rigid body motions in the interior 
domain. We have seen in (32) that this information is lost when trying to prove ellipticity 
ofb(-, ■), but can be reconstructed by adding appropriate terms to b(-, •). We stress that the 
radiation condition (6e) fixes the rigid body motion in the exterior fi ext , see also Section 2.3. 
Since the interior and exterior solution are coupled via equation (28b) this information is 
transferred by (28b) from the exterior to the interior. Thus, adding terms to &(■,■) that 
satisfy (28b) for fixed test-functions seems to be a natural approach. □ 



3.5. Proof of Theorem 2. Let r 1 , . . . , r D be a basis of the rigid body motions TZd and let 
LTo : L 2 (T) — > V®(Eh) denote the 2v 2 -projection. We shall use the observation from the proof 
of Proposition 5 that assumption (33) is equivalent to the fact that rL^r- 7 ), for j = 1, . . . , D 
are linearly independent. 

Proof of Theorem 2 for d = 2. Let 



Mo)< '-(!)• 

denote the canonical basis of TZ2, and let a±, «2, «3 G K fulfill 
(53) ailVr 1 ) + a 2 II (r 2 ) + a 3 n (r 3 ) = 0. 

We stress that n (r 1 ) = r 1 and n (r 2 ) = r 2 . For E G we get 

,3m _ _ 1 (~ j E x 2 dT x 



where s E = (sf, sf) T denotes the midpoint of a boundary element E. Therefore, (53) can 
be written as 

(54) (S) + " 3 ("//)= forall£e ^ 

Altogether, we thus obtain a^s E = ct^s E ' for all E, E' G ££, which can only hold for a 3 = 0. 
This implies o^r 1 + a^v 2 = and hence a± — = Therefore, Ho(r^), j — 1, . . . , 3 = D, 
are linearly independent which is equivalent to (33). □ 

Proof of Theorem 2 for d = 3. Let 

r 1 := |0 I ,r 2 := (l \ , r 3 := | | , r 4 := | Xl I , r 5 - 
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denote the canonical basis of IZ3. We stress that Uq^) = r J for j = 1, 2, 3, and 
n (r- 4 )U = I sA , U (r 5 )\ E = (-si) , n (r 6 )| s = 



*3 



for all faces E G ££, where s B = (sf , sf, sf) T G M 3 denotes the center of mass of an element 
E G E\. The main ingredient for the proof is the geometric observation of Lemma 7 from 
the Appendix: There are at least three elements A,B,C G E\ such that the corresponding 
centers of mass a, b, c do not lie on one line. Let a±, a 2 , 03, 04, as, «6 G K. fulfill 

« 1 n (r 1 ) + a 2 n (r 2 ) + a 3 n (r 3 ) + a 4 n (r- 4 ) + a 5 n (r 5 ) + a 6 n (r 6 ) = 0, 

which is equivalent to 

/aA f-sf 
(55) U + S f - S f 

W V flf 

for all E £ £fr. We take the third equation of (55) for the three elements A, B, C corre- 
sponding to a, b, c and get 






= 0, 

which is only satisfied if 03 = 05 = «6 = or if the vectors 

/l\ /aA /-aA 
(56) I 1 , I 62 I I — b\ are linearly dependent. 




Case 1 (a 3 = a 5 = a e = 0). We insert a 3 , a 5 , a 6 into the first two equations of (55) and 
infer for all E G S\ 

\ + a 4 ( s £? 2 J — or equivalently 



« 4 ^0 = « 4 f s J^ for all £7, £" G £[, 



which can only hold if a\ = 0. Otherwise, all centers of mass would lie on a straight line 
which would contradict the choice of A, B, C. Therefore, we first get 03 = 04 = = = 0, 
and this also implies a\ = a 2 = 0. 

Case 2 ((56) holds). There exist constants 7', 5' G R such that 

(57) U I = 7 ' [ &2 I + 5' 







r «2 




1-1 
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Next, we take the first equation of (55) for the three elements corresponding to a, b, c and 
get 

'\ -a 2 a 3 \ /aA 
1 -b 2 6 3 « 4 = 0, 
1 -c 2 c 3 y \a 6 J 

which is only fulfilled if ot\ = = a§ = or if the vectors 
(58) 



and infer for all E G 




-a 2 \ 






-b 2 


3 


are linearly dependent 


-c 2 







Case 2a («! = a 4 = a 6 = 0). We insert a^a^a^ into equation two and three in (55) 

a2 ) + a 5 ( e 3 J — or equivalently 

« 5 = « 5 ( s s l) f o r a11 # e ^, 

which implies a 5 = as in Case 1. Arguing as above, we first get «i = a 4 = a 5 = a e = 
and finally also a 2 = a 3 = 0. 

Case 2b ((58) holds). There exist constants A,/iGl such that 



a 2 \ 




f«3 


62 


HI 




c 2 y 







Together with (57), we get 
a = a 3 ]^\ J + ' 

where 7 = 7' A and 5 = '-f'fi + 5'. This means that a, 6, c lie on one line which contradicts 
our choice of the elements A,B,C. In particular, case 2b cannot occur. 

Altogether, we have shown a± = a 2 = a 3 = = a$ = a§ = in (55), and therefore 
the orthogonal projections IL^r- 7 ), j = 1, . . . , 6 = D, are linearly independent. Since this is 
equivalent to (33), we conclude the proof. □ 





4. Johnson-Nedelec coupling 

This section deals with the Johnson-Nedelec coupling, see e.g. [JN80, ZKB79] for linear 
Laplace problems and [GHS12, Stel2] for linear elasticity problems. In contrast to [GHS12] 
resp. [Stel2] we avoid the use of interior Dirichlet boundaries resp. an explicit stabilization 
of the coupling equations. The derivation of the variational formulation (59) of the Johnson- 
Nedelec coupling and the proof of equivalence to the model problem (6) are done as for the 
Laplace problem, see e.g. [AFF + 12, GH95] for the derivation. 
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4.1. Variational formulation. The Johnson-Nedelec coupling reads as follows: Find 
(u, 4>)eH = H l {tt) x if- 1/2 (r) such that 

(59a) (2le(it) , e(v)) n - (0 , «}r = (/ , v) n + (</> , v) r 

(59b) , (| - £)u + 5J0)r = , (1 - S)«o)r 

holds for all (v, if)) G "H. Note that the second equation of the Johnson-Nedelec equations 

(59) is the same as for the symmetric coupling (28). We define a mapping i:HxH->K 
and a continuous linear functional F G W by 

(60) b{(u, </>), (v, V)) := (Oe(u) , e{v)) a - (<j> , «) r + <*/> , (f - £)« + 5J0)r 
as well as 

(61) V) := (/ , v) n + (0 O , v) T + (V , (| - £)™o}r 

for all (it, 0), (i>, ?/>) G %. Problem (59) can equivalently be stated as follows: Find (it, <fi) G 
% such that 

(62) 6((«,0),(«,^)) = F(«,^) holds for all (v, t/>) G "H. 
We infer from (60) that 

(63) 6((r,0),(r,0)) =0 for r G 7^. 

Therefore, the mapping b(-, •) cannot be elliptic, and we proceed as in Section 3 to prove 
well-posedness of (59) and its Galerkin discretization. 

4.2. Main result. According to [SW01], there exists a constant 1/2 < c& < 1 such that 

(64) +&)v\\x-i < c A \\v\\ v -i for all v G H 1/2 (r), 

where ||v|||j-i = ( < %3~ 1 v , v) denotes an equivalent norm on H 1/2 (T) induced by the inverse 
of the simple-layer potential. The following theorem is the main result of this section. 

Theorem 6. Let c& < 1 denote the contraction constant (64) of the double-layer potential 
and assume that 2c mon > c$(3A cxt + 2// ext ). Then, the assertions of Theorem 1 hold for the 
Johnson-Nedelec coupling accordingly. 

4.3. Auxiliary results. We stress that the results of Section 3.2-3.3 also apply to the 
Johnson-Nedelec coupling without further modifications. Additionally, the proof needs some 
properties of the boundary integral operators and some results from the works [OS11, Stel2], 
which are stated in the following. First, we introduce the interior Steklov-Poincare operator 
6 : H X/2 {V) -> if^ 1/2 (r) defined by 

6 :=5J- 1 (i + £), 

see e.g. [HW08]. Note that and & are still defined with respect to the exterior Lame 
constants A cxt ,/i ext . We use the estimate 

||(! + #HIIj-i <c&(&w,w) v for all w G H 1/2 (T) 

from [OS11, Stel2], which involves the contraction constant (64) of the double-layer potential 
.ft. The last estimate yields 

(65) (x, (± + R)w)r < + &)w\\fo-i\\x\\fo < V^i&w, w) T \\xh for all (it;, *) G 
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For w G we next introduce the splitting 

(66) w° .— w — w B , 
where w D G H (fi) is the unique weak solution of 

div<x cxt (™ D ) = infi, 

D -p 

w = w on 1 . 

Then, there holds w°\r = as well as the orthogonality relation {<r ext (w D ) , e(w°))n = = 
(<r ext (tu°) , e(w D ))n. Consequently, we see 

(67) (<x cx » , e(w)) u = (<r cxt (™ D ) , e(w D )) n + (<r cxt (™°) , e(u>°)) n . 
Moreover, w D fulfills 7 1 nt, u; D = &w D . Together with Betti's first formula (22), we infer 

(68) (<r cxt (™ D ) , e(w D )) n = (^ti; D , w D ) r = (&w D , w D ) r . 

4.4. Proof of Theorem 6. Note that Proposition 3 holds true with &(■,■) resp. F(-) 
replaced by definition (60) resp. (61). We define the nonlinear operator 33 : T-L — > H* by 

(»(u h ,0 h ), (■,■)) :=6((u fc) (■»■))■ 

Step 1 (Lipschitz continuity o/53). Arguing as in (51) in the proof of Theorem 1, we 
prove Lipschitz continuity of 93, where the Lipschitz constant Cn p > depends only on 21, 
A ext , fi ex \ and Q. 

Step 2 (Strong monotonicity o/Q3). We have to prove that, for all (v,h, 0/J, i v hi V'/J ^ ^> 

(69) {^(u h ,(j> h ) - %(v h ,il> h ) , (u ft - « h ,0 fc - VJ> > C mon ||(u A - w fc ,0 h - 

To abbreviate notation, let (wh, Xh) := ~~ v h, 4>h ~~ tfth)- By use of monotonicity (23) of 
21, we see 

(2*(u h , 0J - Vh) , t»fc>r 

= (2lu h - 2lv ft , w h ) n - ( Xh , w h ) r + (x h , (§ - &)w h + ^Xh)r 

D 

+ ^|(^,(i-^, + 2J X/l ) r | 2 

3=1 

> tw||e(Wh)||£a (n) - (Xfc, (| + ^)wh)r + (Xft, 

D 
3=1 

=: /x - J 2 + J 3 + h. 

Next, we use the splitting (66) for Wh = w° + w D . Together with (67) and (25), where 
d = 6A cxt + 4/i cxt , we get 

/i > , e(w h )) n = -— (<7 ext (™) , e(w u )) n + -—((T ext (w u ) , e(w u )) n 

C-l C-i C-l 

=: iii + ii2- 
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Estimate (65) and Young's inequality yield for 5 > 



With the last inequality and (68), we get 

h< 6 -c^\w»),e(w»)) n + ^\\ X Jl 
Now, we can further estimate the terms I\ — Ii + by 

/i - h + h > hi + (^f " (cr-(« D ) , e(™ D )) n + (l - ^) HXfcllS, 

- (if ~ l c *) ' e ( w *)> n + ( x ~ ir) <x* ' 

where we used (67) again. The assumption 2c mon > cg(3A ext + 2/x ext ) is equivalent to 
Cmon/Cx > c.s/4 with C x = 6A cxt + 4/i ext . Therefore, there exists 5 > such that C := 
min{c mon /Ci — C&S/2, 1 — > 0. Altogether, we infer with (67) and (26) 

D 

h - h + h + h > C<(V xt (™,) , e(w h )) n + ( Xh , Kx h )r + ^ |<£ j , {\ - &)w h + K Xh )r 

> C(||e(^)||^ (n) + , 9Jx h >r + ^ |{^ , (§ - &)w h + K Xh )r\ 2 ) 

= Ci{w hlXh W>CC-^J{w h , Xh )\?u, 

where C = C min{l, C 2 }. The constant C mon := CC^ > depends only on Q, 21, y , and 
on the Lame constants A cxt ,/i ext . 

Step 3 (Unique solvability and Cea lemma). This step is essentially the same as Step 3 
in the proof of Theorem 1. We thus omit the details. 

Remark, (i) In the linear case 21 = cr mt , we may also use an estimate from [Stel2] in 
Step 2 of the proof of Theorem 6 and replace the assumption 2c mon > c$(3A ext + 2/x cxt ) from 
Theorem 6 with 

rj := min{A int /A cxt ,/i int //i cxt } > ^. 

(ii) The assumption 2c mon > c$(3A cxt + 2/i cxt ) ; is an assumption on the monotonicity con- 
stant c mon and the Lame constants A ext ,/i cxt in the exterior domain. As we have seen for 
the symmetric coupling the assumption c mon > suffices to prove unique solvability. Since 
the J ohnson-Nedelec coupling is equivalent to the model problem, we stress that at least 
the continuous formulation of the J ohnson-Nedelec coupling equations is uniquely solvable. 
In [OS11], Of and Steinbach have shown that the J ohnson-Nedelec coupling equations may 
become indefinite (and hence non-elliptic) for special choices of the model parameters. How- 
ever, the numerical experiments from [AFF + 12] show at least numerically that the Laplace 
transmission problem also allows for unique Galerkin solutions in the indefinite regime. 
(Hi) Assume a nonlinear Hencky- Von Mises stress-strain relation, i.e. the operator from (27), 
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with fi(-) > a > and //(•) < Kd/2 — for some a, (3 > 0. T/ien we may replace the as- 
sumption 2c mon > c$(3A ext + 2/x ext ) /rom Theorem 6 with 



where r] := xmn{iaf xeR+ {(K - 2/dfx(x)}/X ex \ inf xeR+ {/i(x)}//i cxt }. □ 



5. Bielak-MacCamy coupling 

In this section we investigate the non-symmetric Bielak-MacCamy one-equation coupling, 
see e.g. [AFF + 12, BM84, CES91] for the Laplace problem. The derivation of the variational 
formulation (70) as well as the proof of equivalence to the model problem (6) essentially 
follow as for the Johnson-Nedelec coupling resp. symmetric coupling, cf. e.g. [AFF + 12, 
CFS97, GH95]. 

5.1. Variational formulation. The variational formulation of the Bielak-MacCamy cou- 
pling reads as follows: Find (it, (f>) G % = H 1 ^) x H~ 1/2 (T) such that 

(70a) (2le(w) , e(«)) n + ((± - , v) T = (f , v) n + (0 O , v) T 

(70b) (i/>, QJ0- w ) r = -(,/,, n ) r 

holds for all (v, iji) G We sum up the left-hand side and the right-hand side of (70) and 
define the mapping b : Ti x H — > K as well as the linear functional F G by 

(71) 6((u, 0), (t>, V)) := (2te(u) , e(«)) n + ((± - , t>) r + , ^ ~ ™}r 
as well as 

(72) F{v, V) := (/ , v) n + (0 O , w) r - (V , n )r 

for all (it, 0), (v,if>) G Then, problem (70) can equivalently be stated as follows: Find 
(it, (p) G H such that 

(73) b((u,<p),(v,ip)) = F(v,^) holds for all (i;, V) G 

As for the other coupling formulations &(•,•) is not uniformly elliptic, and unique solvability 
cannot be shown directly. We follow the ideas of Section 3 resp. Section 4 to overcome these 
difficulties. Moreover, with 6jn(-, •) denoting the mapping defined in (60), we stress that 

(74) 6 JN ((u, </>), (u, 0)) = b((u, 0), (u, 0)) for all (ti, <f>) G ft. 

Thus, there is a strong relation between the one-equation Bielak-MacCamy and Johnson- 
Nedelec coupling. In fact, for linear and symmetric 21 : L 2 (Q) — > L 2 (Q), there holds 

b m ((u, (f>), (v, VO) = 6(( v > V>)> 0> 0)) for a11 K 0), K e W. 



5.2. Main result. As in Section 3.2, we add appropriate terms to &(•, •) to tackle the rigid 
body motions in the interior domain Vt. In particular, we use (70b) to stabilize the linear 
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form b(-, ■). We stress that Proposition 3 holds with b(-, •) resp. F(-) replaced by 

D 

(75) b((u, 0), (v, 0)) := b((u, 0), (v, */>)) + , 2J0 - u) r (£ , - </>) r , 

(76) F(v, <#) := F(«, V) - ' u °> r ^ ' " 

for all (it, 0), (v, 0) G Furthermore, the assertions of Proposition 5 also hold true if (46) 
is replaced by 

(77) <£,9J0-u)r for (u, 0) G ft. 

With these observations, Theorem 6 holds true for the Bielak-MacCamy coupling. Details 
are left to the reader. 

Remark. Our techniques developed in Section 3-5 can also be used for the (quasi-) sym- 
metric Bielak-MacCamy coupling schemes [BM84, GHS12] applied to nonlinear elasticity 
problems. □ 

Appendix A. Elementary geometric observation 

Lemma 7. Let d = 3 and E\ be a regular triangulation of the closed boundary T = dQ into 
flat surface triangles. Then, there are at least three triangles A,B,C G E\ such that the 
centers of mass a, 6, c corresponding to these elements do not lie on one line, i.e. c — a ^ 
{t(b-a) : teffi}. 

Proof. We argue by contradiction. Assume that all centers of mass lie on one line g. Let 
x G /C£ denote an arbitrary node of the triangulation E\. Recall that V = dQ is the 
closed boundary of the polyhedral Lipschitz domain Q. Therefore, there are k > 3 triangles 
Di, . . . , Dk G £fr such that a; is a corner of Dj for j — 1, . . . , k. Moreover, let {x, x l , . . . , x k } 
denote the set of all nodes of the triangles Dx, . . . , D^. We stress that we can permute the 
indices of x l , . . . , x k and the indices of Dx, . . . , such that Dj = conv {x, x\ x^ +1 }, where 
we define x k+1 := x 1 and x k+2 := x 2 , see Figure 1 for an illustration in the case k = 6. Let 
s 1 , . . . , s k denote the centers of mass of the triangles Di, . . . , Dk- Because of our assumption 
that all centers of mass lie on one line g, we infer that 

, , , , x j+2 + x j+l + x x j+l + x j + x x j+2 - x j 
s J — S J = = 

3 3 3 

is proportional to the directional vector d ^ of the line g. Therefore, a? J+2 — x^ is also 
proportional to d. Moreover, we observe td = s J+3 — s J+2 + s J+1 — s J = (x^ + — x^)/3 for 
some tel. By iterating this process, we get with appropriate t m ,t n G K 

2m 

t rn d = ^(-1)V = {x 2m+1 - x 1 )/^ for all m with 2m < k, 
j'=i 

2n-l 

t n d = (-l) j+1 s j = {x 2n - x 2 ) /3 for all n with 2n - 1 < Jfe. 

J'=2 
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Figure 1. For any node x G 1C\ in a regular triangulation E\ of a closed 
boundary T, there exist k > 3 different nodes x 1 , . . . , x k G /C^ and triangles 
D\, . . . , Dk such that Dj = conv {cc, a? J , a?- 7-1-1 } for all j = 1, . . . , k. Moreover, 
there holds Dj Pl-Dj+i = conv {x, x^ +1 } for j = 1, . . . , k with cc fc+1 = a? 1 . Here, 
an example for k = 6 is shown. 

Altogether, we see that all nodes with even indices lie on one line h, and all nodes with odd 
indices lie on one parallel line f, i.e. 

x 2j e {x 2 + td : t e R} =: t) and 
x 2 ^ 1 G {x 1 +td : t G M} =: f 

for all j G N with 2j < k resp. 2j — 1 < k. For the remainder of the proof, we distinguish 
whether k is odd or even. 

Case 1 (k is odd). The observations above show that x 1 , x 3 , . . . , x k G f and x 2 — x k = 
x k+2 — x k = td for some t G R. Then, x 2 6 f and since f and t) are parallel, there holds 
h = f, which means that all nodes x l , . . . , x k lie on one line. This contradicts a regular 
triangulation. 

Case 2 (k is even). If h = f we can argue as in Case 1. Otherwise f) ^ f, and we stress 
that all edges x^x^ +1 of the triangles are connected, i.e. 

{{x j } ifj = n+l 
{x j+1 } ifj + l = n 
otherwise. 

Moreover, every edge x^x^ +1 connects the lines t) and f. Thus, we can infer that there are 
two edges x m x m+1 , x n x n+1 which intersect each other, i.e. 

x m x m+i n x n x ii+\ = j^} w ^ y ^ x n and y ^ x n+1 , 
see also Figure 2 for an illustration. Altogether this contradicts a regular triangulation. □ 
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Figure 2. Let x e /C£ be an arbitrary node in the regular triangulation S,\ 
and let x 1 , . . . , x k denote the neighboring nodes. Under the assumption that 
all centers of mass of the triangles in E\ lie on a line g, the proof of Lemma 7 
unveils that the nodes £C 2j_1 lie on a line f and the nodes x 2 i lie on a line 
f), which is parallel to f. Since all a?- 7 are connected by the segments x^x^ +1 , 
we can conclude that there are indices n, m with 2n < k, 2m < k such that 
x 2n x 2n+1 H x 2m x 2m+1 = {y} and y £ K\. 
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